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Abstract 

We study compact stars that contain quark matter. We look at the effect 
of color superconductivity in the quark matter on the nuclear-quark matter 
transition density, nrass-radius relationship, and the density discontinuity 
at the boundary between nuclear and quark matter. We find that color 
superconducting quark matter will occur in compact stars at values of the 
bag constant where ordinary quark matter would not be allowed. We are 
able to construct “hybrid” stars with a color superconducting quark matter 
interior and nuclear matter surface, with masses in the range 1.3-1.6 Mg and 
radii 8-11 km. Our results are consistent with recent mass-radius limits based 
on absorption lines from EXO0748-676. 


1 Introduction 


If matter is compressed far enough beyond nuclear density then there is a transition 
from nuclear matter to quark matter. It is becoming widely accepted that quark 
matter will typically be in a color-superconducting phase I | i I- 0- in which 
the quarks near the Fermi surface form Cooper pairs which condense, breaking the 
color gauge symmetry (for reviews see Ref. 0). The pairing pattern favored at 
sufficiently high density is the color-flavor locked (CFL) phase in which up-down, 
down-strange, and up-strange Cooper pairs all form, allowing quarks of all three 
colors and all three flavors to pair |fj. 

One of the most likely locations for quark matter in nature is the interior of 
compact stars, where pressure due to gravity drives the density above nuclear 
density, and the temperature is low compared with nuclear/quark energy scales. 
Various possible signatures of color superconductivity in compact stars have been 
studied, mostly focusing on transport properties (for a recent review, see Ref. [[!]]). 
Although the effect of unpaired quark matter on the compact star mass-radius 
relationship is an active area of research I I 0, the consequences of color 


superconducting quark matter have not yet been investigated. 

The contribution of color superconductivity with gap A to the pressure p of quark 
matter is of order /r 2 A 2 , which is dominated by the leading order /i 4 contribution 
from the Fermi sea. However, quark matter must also pay a free-energy cost, the 
bag constant R, relative to the confined vacuum: 
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If the bag constant is large enough so that nuclear matter and quark matter 
have comparable pressures at some density that occurs in compact stars, then the 
superconducting gap A may have a large effect on the equation of state and hence 
on the mass-radius relationship of a compact star. A similar observation has been 
used recently to show that the region of model space where strange quark matter is 
absolutely stable is influenced by superconductivity |Ll . 


In this paper we study the gross structure of compact stars, taking into account 
the possibility of a CFL quark-paired phase of quark matter, as well as unpaired 
quark matter (UQM). We treat the quark matter as a Fermi sea of free quarks with 
an additional contribution to the pressure from the formation of the CFL condensate. 
We treat the nuclear matter as consisting only of protons, neutrons, and electrons, 
and either obtain the nuclear equation of state from the Walecka model of nuclear 
interactions, or use the APR98 equation of state 


12 


In order to fold the effects of color superconductivity, we allow all the phases to 
compete with each other, selecting the highest pressure phase at each value of the 
quark chemical potential. If local electrical neutrality is imposed, this leads to sharp 
interfaces between the different phases. If a globally neutral inter-penetration of two 
charged phases is allowed, this leads to mixed phases. We study both possibilities. 
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We survey the parameter space of the quark matter, covering a range of values of 
the bag constant B, the strange quark mass rn s , and the color superconducting gap 
A. In section [2] we explain our calculation of the equations of state for nuclear 
and quark matter. In section we obtain the resultant mass-radius relationships, 
paying particular attention to the maximum masses that can be obtained. Section 
f| presents our conclusions. 


2 Equations of state 


For nuclear matter we use the Walecka equation of state, which allows us to 
calculate the pressure for any quark chemical potential /i and electron chemical 


potential /z e , and hence to construct mixed phases. In section |3.3| we will also show 
results for the APR98 equation of state, which is obtained using non-relativistic 
variational methods starting from a Hamiltonian that reproduces known nucleon- 
nucleon scattering data. For nuclear matter at very low densities we use the 


tabulated Negelc-Vautherin |13| and Baym-Pethick-Sutherland JT4| equations of 
state. 


2.1 Walecka equation of state 


We use the Walecka model as described in in and calibrated in ]16|. The free 
energy density is 
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where 


e n {k) = \jk 2 + m* N 2 + g uN u - \g pN p , (2.2) 

e p (k) = \jk 2 + m* N 2 + g uN u + \g pN p , ( 2 . 3 ) 

are the neutron and proton single particle energies in the mean held approximation. 
The corresponding Fermi momenta kp n and kp v , which minimize the free energy at 
fixed baryon and electron chemical potentials, are given by solving 

£n(kpn) Hn > 

£p(kp p ) Up , 

where weak interaction equilibrium sets /i p = ii n — /i e , and 

m 2 a a = g aN (n s (k Fn ) + n s (k Fp )) - ^ , 
mlcu = g uN (n(k Fn ) + n(k Fp )) , 
m 2 p = \g p (n(k Fp ) - n(k Fn )) . 


(2.4) 
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The nucleon number density n and scalar density n s for nucleons with Fermi 
momentum kp are 


where 
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is the nucleon effective mass, which is reduced compared to the free nucleon mass m n 
due to the scalar field a, taken to have m a = 600 MeV. The scalar self-interaction 
term is 


U(a) 


^ m N {g cN af + ^{g aN a) 4 


( 2 . 8 ) 


where b and c are dimensionless coupling constants. The five coupling constants, 
g^Ni g^N, g P N, b , and c, are chosen as in Ref. |L6| to reproduce five empirical 


properties of nuclear matter at saturation density: the saturation density itself is 
no = 0.16 fm -3 ; the binding energy per nucleon is 16 MeV; the nuclear compression 
modulus is 240 MeV; the nucleon effective mass at saturation density is 0.78?7 Ttv; 
and the symmetry energy is 32.5 MeV. 

The charge density in nuclear matter is 


rn _ '^nuclear (r) ^ 

nuclear 0 V ^ ^ / 

Ofl e 

which is just the number density of protons minus that of electrons. In bulk matter 
one requires Q nuclear = 0, but not in mixed phases (Section pT5|) . 


2.2 Unpaired quark matter (UQM) equation of state 

In noninteracting unpaired quark matter, and neglecting the light quark masses, the 
free energy density is 
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where the Fermi momenta are 


”1 = Pi ~ m l where Pu = p-\p e , 


u d ~ Pd ~ m d where p d = p + \p 


u s = P-l ~ m s where p s = p 
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Differentiating with respect to p e , we obtain the charge density 
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2.3 CFL quark matter equation of state 

We describe the CFL phase using the free energy 

^CFl(/L Pe) = ^CfTXu) + ^CFl(/L Pe) + ^ 
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The contribution to (|2.13| ) from the quarks is || [17 
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where the quark number densities are n u = n d = n s = (z/ 3 + 2A 2 p)/7r 2 and the 
common Fermi momentum is 


v = 2/z- \/p 2 + 
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(2.15) 


The first two terms give the free energy of the noninteracting quarks, while the third 
term is the lowest order (in powers of A/p) contribution from the formation of the 
CFL condensate. 

The contribution to (|2.13|) from the Goldstone bosons arising due to breaking 
of chiral symmetry in the CFL phase is denoted 00p L (/i,/r e ). The effective theory 
describing the octet of mesons has been studied extensively in earlier works JT8| . 
When the electron chemical potential exceeds the mass of the lightest negatively 
charged meson, which in the CFL phase is the n~, these mesons condense |TP], |20fl . 
The free energy in this case is given by 
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where the parameters are JT8| 
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Finally, the contribution to ( |2.18| ) from electrons is 
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Unlike the UQM phase, in the CFL phase there is a gap in the quark excitation 
spectrum, and the lightest charged excitations correspond to pions and kaons. The 
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charge susceptibility in this phase is determined by the effective theory for these 
collective modes. The electric charge density (carried by the pion condensate) 
induced by a electron chemical particle is given by [i2Q| 
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Meson condensation can occur in the CFL phase even in the absence of an 
electric charge chemical potential. Bedaque and Schaefer |H| have shown that the 
strange quark mass introduces a stress on the CFL state and that might result in the 
condensation of K° mesons. Condensation occurs when m 2 s /2n > rn K o, where m^o 
is the mass of the K° meson in the CFL phase. The free energy due to A'°-condensed 
phase is 
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where the kaon mass 18 
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From Eq. ( 12.20 ) we see that the free energy due to A'°-condensation is an order 
m 4 effect and thereby small compared to the A 2 /i 2 contribution to the free energy 
for A ~ 100 MeV. For this reason we neglect K° condensation in this study. 


2.4 Is color superconductivity important for bulk structure? 

We see from ( p,14| ) that color superconductivity contributes 0(yU 2 A 2 ) to the free 
energy, which is small relative to the kinetic energy density which is 0(/i 4 ). This 
well-known suppression is a consequence of the fact that pairing is a Fermi surface 
phenomena and the superconducting gap is usually small compared to the chemical 
potential. Naively, this would lead us to conclude to that superconductivity will not 
greatly affect the equation of state of quark matter. If this were true, we should 
expect that the mass-radius relation of neutron stars containing superconducting 
quark matter would be nearly identical to those constructed in earlier works 
wherein the role of superconductivity was neglected. However, the situation is more 
complicated. In the bag model description of quark matter, the free energy gets 
an additional contribution due to the bag constant. The kinetic pressure and bag 
pressure cancel when the quark chemical potential has value 
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Thus, for a given B, there is a narrow window in quark chemical potential in 
which the pairing contribution to the pressure is dominant. In the vicinity of /i 0 , 
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Bag constant 

CFL 

chemical 

transition 

nuclear 

CFL 

B 1 / 4 

B 

gap 

potential 

pressure 

density 

density 

(MeV) 

(MeV/fm 3 ) 

(MeV) 

(MeV) 

(MeV/fm 3 ) 

(?Lat) 

(^sat) 

190 

169.6 

0 

422.5 

111.5 

3.464 

5.833 

190 

169.6 

50 

408.2 

88.68 

3.158 

5.405 

190 

169.6 

100 

365.4 

33.79 

2.149 

4.295 

170 

108.7 

0 

352.6 

21.62 

1.805 

3.297 

170 

108.7 

50 

338.4 

10.77 

1.382 

3.033 


Table 1: Properties of the nuclear-quark phase transition for various bag constants 
B and color-superconducting gaps. Nuclear matter is treated using the Walecka 
model. The size of the gap has a significant effect on the pressure at which the 
phase transition occurs and the densities of the two phases there. 


superconductivity will therefore make a significant contribution to the equation of 
state of quark matter. For B 1 ^ in the range 150-200 MeV (B = 66-210 MeV/fm 3 ), 
we find /x 0 — 320 — 400 MeV. This is an interesting range of chemical potentials 
because the phase transition from nuclear matter to quark matter typically occurs 
here. Further, and perhaps more importantly, the pairing contribution to the 
pressure of CFL quark matter, Pa = 3A 2 / u 2 /7t 2 , can be comparable or larger 
than the pressure in the nuclear phase at the same baryon chemical potential. 
Superconductivity will thereby significantly influence the critical chemical potential 
at which the transition from nuclear to quark matter occurs. This is clear from 
Fig- 0, where the transition from NM to QM occurs at a much lower pressure for 
CFL QM (a— >b) than for unpaired QM (c— >d). The critical chemical potential for 
the transition for different values of the bag constant and the superconducting gap 
are shown in table 0 


2.5 Mixed phases 


In the preceding discussion we have enforced local charge neutrality in the nuclear 
and quark phases. We have neglected the possibility of having a mixed phase, at 
finite fi e , containing positively charged nuclear matter co-existing with negatively 
charged CFL matter pH] . Such a possibility was considered in Ref. |22|] where the 
bulk free energy difference between the homogeneous phases and the heterogeneous 
mixed phase was calculated. The free energy difference between the homogeneous 
phase and the mixed phase was found to be quite small. If one accounted for the 
additional surface energy cost in the mixed phase, it was found that even modest 
values of the surface tension, onm-cfl ~ 30 MeV/fm 2 were sufficient to disfavor the 
mixed phase. Nonetheless, we consider this possibility in this work for two reasons. 
First, the surface tension between the nuclear matter and CFL matter is poorly 
known. Second, in the case of unpaired quark matter, allowing for a mixed phase 
has been shown to significantly affect the equation of state over a wide range of 
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pressures and consequently modify the mass-radius relation. This indicates that the 
superconducting case warrants investigation. 

The procedure to construct the mixed phase between nuclear matter and CFL 
matter was outlined in Ref. |22|. We follow the same prescription here but briefly 
note some salient features which distinguish the nuclear-CFL mixed phase from 
the mixed phase between nuclear matter and unpaired quark matter. The volume 
fraction of the nuclear and quark phases in the mixed phase is determined by the 
requirement of global charge neutrality. Denoting the charge density of the nuclear 
phase by Q nuclear and the charge density of the CFL phase as Qcfl, the volume 
fraction of the CFL phase is 


X = 


Q 


nuclear 


VcFL 

C Qnuclear Q CFL 


(2.23) 


The energy density of the mixed phase is the volume-weighted average of the 
individual energy densities and is given by e = x 6 cfl + (1 — x)e nU ciear- From 
Eq. ( 2.23|) , we see that, if at fixed electron chemical potential the negative charge 
density of the quark phase is large its volume fraction will be correspondingly 
smaller. This is important because the charge susceptibilities of the normal and 
superconducting phases are quite distinct. In the normal phase it is easy to furnish 
electric charge since there is no gap in the spectrum for quarks. The electric charge 
densities of the nuclear matter (Walecka), unpaired quark matter, and CFL quark 
matter phases are given in Eqs. (ED> O, O- 

For typical quark and electron chemical potentials encountered in the neutron 
star context one finds that the charge density in the normal phase is significantly 
larger. For example, when /i = 400 MeV, m s = 150 MeV, /j, e = 100 MeV and 
rriK = 30 MeV, the charge density in the normal phase Quqm = —0.32 frn -3 , while 
in the CFL phase Qcfl = —0.09 frn -3 (note that neutral CFL QM has /i e = 0 
17| , so the charge density in CFL is due to the density induced by /i e ). From 


the preceding arguments this implies that the volume fraction of the CFL phase 
in the mixed phase will be significantly larger. For the same reason, although the 
pairing contribution to the free energy itself is small, the equation of state of the 
CFL-nuclear mixed phase is considerably softer. 

The equation of state for the CFL-nuclear and UQM-nuclear mixed phases are 
shown in Fig. [I]. We see that if sharp transitions occur, then the occurrence of 
color superconductivity leads to higher energy density at low pressure because 
the NM—>CFL transition (a —> b) occurs at lower pressure than NM—^unpaired 
(c—>d), but lower energy density at high pressure, because CFL has a lower energy 
density than unpaired QM. On the other hand, if mixed phases occur, then color 
superconductivity leads to a higher energy density up to high pressures, since the 
NM+CFL mixed phase {A—> B) and CFL phase that follows it both have higher 
energy density than the NM+unpaired mixed phase (C —> D) over a wide range of 
pressures. 
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Figure 1: The equation of state for Walecka nuclear matter and quark matter with 
T? 1 / 4 = 185 MeV (B = 153 MeV/fm 3 ) and m s = 200 MeV. We show normal unpaired 
quark matter (thin lines) and color flavor locked superconducting quark matter with 
A = 100 MeV (thick lines). We show neutral phases with sharp transitions (solid 
curves) and mixed phases (dashed curves). Differences between the Walecka model 
equation of state and the APR98 equation of state can be inferred from the figure 
inset. 

A mixed phase involving CFL and UQM is also allowed in principle. Such a 
mixed phase is characterized by a negatively-charged CFL phase coexisting with a 
positively-charged UQM phase. We find that such a mixed phase has lower free 
energy (even in the absence of surface and Coulomb effects) only if the negatively 
charged mesons (kaons) condense in the CFL phase. This requires that /i e > mx- 
On the other hand, the UQM phase is positively charged only when /i e < m 2 / 4/i. 
Further, the difference in free energy between this phase and the pure CFL phase 
is only of order m 4 and it does not greatly affect the equation of state. For these 
reasons we do not consider such a mixed phase in our study of the structure of 
the compact star. Nonetheless, the existence of such a mixed phase might have 
important consequences for transport and cooling phenomena especially in pure 
quark stars. 
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3 Compact Star Structure 


To determine the mass and radius of the compact object for a given value of the 
central pressure, we must solve the Tolman-Oppenhcimer-Volkov (TOY) equations 


dP 

dr 

dM[r) 

dr 


-G M{r) e(P) f P\ 

r 2 C 2 V e > 
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Anr 3 P \ 
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where P = P(r ) and the equation of state specifies e(P), i.e., the energy density as a 
function of the pressure, and M{r ) is the total energy enclosed within radius r. For 
a given central pressure, P(r = 0), the above equations can be easily integrated out 
to the surface of the star, where P = 0, to obtain the mass and radius of the object. 
By varying the central pressure it is possible to obtain the mass radius relation 
predicted for a given model description of the matter equation of state. The focus 
of this section is to employ the equations of state described in the previous section 
and deduce the corresponding mass-radius relationship. 


3.1 Ad-R relationship for uniform stars 


First, we review the mass-radius relationship for the simple cases of stars made of 
pure nuclear matter or pure quark matter. The two solid curves in Fig. ^ are for 
pure nuclear matter stars. They follow from solving the TOV equation at high 
bag constant, so quark matter is highly disfavored and does not occur. There is 
one for each of the two nuclear equations of state that we use in this paper: the 
Walecka model (Section |2.1|) and the APR98 tabulation [T2]|. They are roughly 
similar, showing the characteristic sharp rise in mass at R ~ 10-15 km as the 
NM equation of state hardens around nuclear saturation density. They also show 
similar maximum masses of about 2 M Q , at central densities of 5 to 10 times nuclear 
saturation density. Despite these similarities there are, as can be seen in the inset 
of Fig. m, quantitative differences. The APR98 equation of state is on average softer 
at low density and stiffer at high density compared to the mean held model. For a 
detailed discussion of how different nuclear equations of state affect the mass-radius 
relationship of neutron stars see Ref. ||25|| . 

For comparison, Fig. also shows (dashed line) the M-R curve for quark matter 
with no pairing (A = 0) at a low bag constant (P 1//4 = 145 MeV, B = 58 MeV/fm 3 ) 
where three-flavor quark matter is favored over nuclear matter all the way down to 
zero pressure, but two-flavor quark matter is less favorable than nuclear matter at 
low pressure. The maximum mass is very sensitive to the bag constant, so the fact 
that this curve also shows a maximum at around 2 M 0 is coincidental. (The mass 
and radius scale as 1 /\/~B |T5|.) Not all the values of M and R that lie on the curves 
are stable. The family of stable configurations is generated by increasing the central 
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Figure 2: Mass-radius relationships for pure nuclear matter (NM) with Walecka 
and APR98 equations of state, and for unpaired quark matter with B 1//4 = 145 MeV. 
The straight dashed line indicates the constraint obtained by recent measurements 

m 


of the redshift of spectral lines from EXO0748-676 


This constraint requires 


model equations of state to intersect the line M/R m 0.15 (M©/km) 


pressure, and obtaining an increasing mass. As soon as the maximum of M(p cen trai) 
is attained, further increases in Pcentrab apparently yielding lighter stars, will in fact 
move onto the unstable branch. This means that the the parts of the Ad-R curves to 
the left of the maxima in Figs. [|, || correspond to stars that are hydrodynamically 
unstable to collapse to a black hole. 


3.2 M-R relationship for hybrid and color-superconducting stars 


The main purpose of this paper is to explore the effect of quark pairing on the 
M-R relationship at values of the bag constant that are consistent with nuclear 
phenomenology. Fig. |3] shows the mass-radius curve for a plausible model of dense 
matter: the Walecka nuclear equation of state, and quark matter with physically 
reasonable values of the bag constant B 1//4 = 180 MeV (B = 137 MeV/fm 3 ) and 


strange quark mass m s = 200 MeV 0. Curves for unpaired (A = 0) and color¬ 
superconducting (A = 100 MeV) quark matter are shown. At these values the 
stars are typically “hybrid”, containing both quark matter and nuclear matter. The 
solid lines in Fig. [3] correspond to stars that either have no QM at all, or a sharp 
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transition between NM and QM: the core is made of quark matter, which is the 
favored phase at high pressure, and at some radius there is a transition to nuclear 
matter, which is favored at low pressure. The transition pressure is sensitive to A, 
for reasons discussed earlier. The dashed lines are for stars that contain a mixed 
NM-QM phase. In all cases we see that light, large stars consist entirely of nuclear 
matter. When the star becomes heavy enough, the central pressure rises to a level 
where QM, either in a mixed phase or in its pure form, occurs in the core. As can 
be seen from the figure the transition density is very sensitive to A. 
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Figure 3: Mass-radius relationships at fixed bag constant B 1 ^ = 180 MeV and 
m s = 200 MeV, for unpaired (A = 0) and color-superconducting (A = 100 MeV) 
quark matter. The mixed phase (dashed) and the sharp interface curves are shown. 
The line labeled “Cottam et al .” indicates the constraint obtained by recent 
measurements of the redshift on three spectral lines from EXO0748-676 |24|. The 


dots labeled p 0 and 2 p 0 on the nuclear matter mass-radius curve indicate that the 
central density at these locations correspond to nuclear and twice nuclear saturation 
density respectively. The top panel uses the Walecka equation of state for nuclear 
matter, and the lower panel uses APR98 (in which case we only consider the sharp- 
interface scenario). 
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The profiles of the maximum mass superconducting stars for different values of 
the bag constant, A = 100 MeV and rn s = 200 MeV are shown in Fig. f|. For 
B l / A = 185 MeV results for the sharp interface (denoted as (s)) and the mixed 
phase (denoted as (m)) scenario are shown. Here the maximum masses correspond 
to 1.33 M q and 1.35 M 0 , respectively. The maximum mass for B 1 / 4 = 175 MeV and 
H 1 / 4 = 170 MeV are M max = 1.44 M 0 and M max = 1.52 M 0 , respectively. Fig. f] 
shows that the typical density discontinuity in the sharp interface scenario is ~ 3 p a . 
It also shows that for smaller values of B , the NM ^ QM phase transition occurs 
very close to the surface of the star (at lower density as discussed earlier). The 
denser exterior regions of these stars (despite a less dense inner core) are primarily 
responsible for the increase in the maximum mass observed as one decreases B. 



Figure 4: Profile of the maximum mass star for bag constant B 1//4 = 
185, 175, 170 MeV with m s = 200 MeV and A = 100 MeV. The mixed phase 
(dashed) and the sharp interface curves are shown. The Walecka model was used 
to describe the nuclear part of the equation of state. 

Fig. [| indicates that in the mixed phase scenario there are no discontinuities in 
the density profile of the star. However, this is not true in general. It is interesting 
to note that even when mixed phases are allowed, there can still be discontinuities 
in energy density within them. In a small range of parameters, we find stars that 
have a crust of nuclear matter surrounding a mixed NM-QM core, but the mixed 
phase has an outer part which is a mixture of unpaired QM with NM, and an inner 
part that is a mixture of CFL QM with NM. At the interface between the two there 
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is a density discontinuity within the mixed phaseQ. 

The results shown in Figs. |3].[| indicate several generic trends: 

1. The stability of stars containing pure QM but without a QM-NM mixed phase 
(i.e., a sharp-interface scenario) depends on the transition density. Color 
superconductivity has a favorable effect on the stability of these objects since 
it lowers the transition density and stiffens the equation of state relative to 
UQM. 

2. In general, mixed phase stars are more likely to be stable since the equation 
of state (or adiabatic index) changes smoothly in this case. For example, at 
A = 0 the uniform phase stars become unstable as soon as quark matter is 
introduced (the A = 0 line slopes down to the left, showing a decreasing mass 
as the central pressure rises) whereas the mixed phase (dashed line) gives a 
stable branch leading to a maximum mass M max ps 1.5 M 0 . 

3. For large A there is very little difference between the sharp interface and mixed 
phase mass-radius curves. This is because the volume fraction of the CFL 
grows very rapidly within the nuclear-CFl mixed phase as discussed earlier 
in section p3| . Consequently, the extent of the mixed phase is reduced and 
the equation of state with in this region more closely resembles the pure CFL 
equation of state as is evident from Fig. [[]. 

4. Although there are visible differences between the mass-radius curves of 
the stars with UQM or CFL matter depending on whether one employs 
the Walecka or the APR98 equation of state, the maximum mass and 
corresponding radius of such a star is fairly independent of the nuclear equation 
of state. This is because the nuclear phase contributes very little to the total 
mass of these stars. These features of the mass-radius relation will be more 
comprehensively studied in section |3^ . 

5. If the transition to superconducting quark matter is constrained to occur at 
or above nuclear density the maximum mass of these stars is ~ 1.4 M 0 . It 
is possible to obtain larger masses (~ 1.6 M 0 ) if the transition is allowed to 
occur at lower density. We elaborate further on this in the subsequent section. 

3.3 Color superconductivity and the maximum mass 

In bag model treatments of quark matter, the bag constant, strange quark mass, 
and the color superconducting gap are unknown parameters. In this paper we 
take a reasonable range of values for B and m s , and study the dependence on 
A of observable features of compact stars such as their mass and radius. The 
resultant predictions can then be used to constrain the CFL color superconducting 
gap. In Fig. ^|we show how the maximum star mass (obtained by varying the central 

1 The inclusion of additional phases, such as the two flavor superconducting (2SC) phase, might 
modify this conclusion. 
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pressure) depends on color superconducting gap A for two different bag constants 
and two different strange quark masses. 

For each value of the bag constant and strange quark mass there are two or three 
curves of M max vs A. The solid curve is for stars with some quark matter, paired 
or unpaired in them. The prominent dot on the curve separates the pure quark 
stars (to the right, at higher gap) from the stars with a QM core and a NM crust or 
mantle (to the left, at lower gap). The dotted curve indicates the heaviest pure NM 
star. This depends on the gap because A affects the point in the M-R plot at which 
quark matter appears. For example, for the equations of state used in the top panel 
of Fig. the maximum NM mass would be 0.8 M 0 , (radius 13.8 km), since that is 
where QM first appears, and the star ceases to be pure NM. For the Walecka model 
that we used, the maximum possible NM mass is about 2.04 M 0 . For the APR98 
equation of state, the maximum possible NM mass is about 2.20 M 0 . In general, a 
large gap favors quark matter, causing the heaviest pure NM star to become lighter. 

(1) Large bag constant, B 1//4 = 185 MeV, B = 153 MeV/fm 3 
(Upper two panels in Fig. ^j). 

(a) Pure NM stars still occur, but their maximum mass drops as the gap grows 
(dotted lines in Fig. ||) as the unstable QM branch cuts off the NM branch. 

(b) Stars with a QM core and NM surface separated by a sharp interface (solid lines 
to the left of dots) become stable at gap A ~ 50-120 MeV, depending on m s . These 
stars have mass M < 1.5 M 0 . 

(c) Stars with a NM-QM mixed phase core and NM surface (dashed lines) occur at 
lower values of the gap, and masses up to around 1.7 M 0 are possible if the strange 
quark is heavy enough. 

(d) Pure QM stars (solid lines to the right of dots) occur at very large gaps. 

(2) Small bag constant, T? 1 / 4 = 165 MeV, B = 96 MeV/fm 3 
(Lower two panels in Fig. |5j). 

(a) Pure NM stars are again cut off by the QM branch. 

(b) Stars with a QM core and NM surface separated by a sharp interface (solid 
lines to the left of dots) have masses up to about 1.6 M 0 . For light strange quarks 
m s ~ 150 MeV, color superconductivity increases the maximum mass attained by 
the stars, but from a lower starting point at A = 0. 

(c) Stars with a NM-QM mixed phase core and NM surface (dashed lines) also have 
masses up to around 1.6 M 0 if the strange quark is heavy enough. 

(d) All QM-containing stars with mass greater than about 1.6 M 0 are pure quark 
stars (solid lines to the right of dots). 

Our overall conclusion is that in this range of values of the bag constant, turning 
on color superconductivity allows hybrid stars to occur with masses up to around 
1.6 Mq. In Section f| we will discuss in more detail whether this is compatible with 
recent observational data, and how it should be used to interpret future observational 
data. 
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Figure 5: The maximum star mass (in solar masses) attained as a function of color 
superconducting gap A. The upper plots are for T? 1 / 4 = 185 MeV, the lower plots 
for H 1 / 4 = 165 MeV. The left plots are for nuclear matter described by the Walecka 
model (in which case a mixed phase can be constructed). The right plots are for the 
APR98 nuclear equation of state for which we have only constructed locally neutral 
phases. Curves for strange quark mass m s = 150 and 300 MeV are shown. The 
solid lines give the radius of stars with QM only (to the right of the dot), or a QM 
core surrounded by NM (to the left of the dot). The dotted lines show the heaviest 
pure NM star that occurs at the given gap A. The dashed lines are for stars that 
include a mixed phase. 
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Figure 6: The radius of the star (in km) when it attains its maximum mass as 
plotted in Fig. [5|. The upper plots are for B 1 ^ = 185 MeV, the lower plots for 
F? 1 / 4 = 165 MeV. The left plots are for nuclear matter described by the Walecka 
model (in which case a mixed phase can be constructed). The right plots are for 
the APR98 nuclear equation of state. Curves for strange quark mass m s = 150 and 
300 MeV are shown. The solid lines give the radius of stars with QM only (to the 
right of the dot), or a QM core surrounded by NM (to the left of the dot). The 
dotted lines show the heaviest pure NM star that occurs at the given gap A. The 
dashed lines are for stars that include a mixed phase. 
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3.4 Color superconductivity and compactness 

By comparing the A = 100 MeV curves of Fig. |3| with the UQM curve in Fig. |2j. 
we see that superconductivity only has a moderate effect on the radius of the star. 
In Fig. we investigate this issue more comprehensively, by plotting the radii of 
the stars whose masses appeared in Fig. [5], ie the radius of the heaviest star at each 
value of A for each equation of state. This confirms that color superconductivity 
does not affect the radii very strongly, but at low to medium gap A < 100 MeV 
it tends to reduce them. This means that increasing A is not simply equivalent to 
decreasing the bag constant, in which case the mass and radius grow together as 
1 /\[B. By comparing Fig. [5] with Fig. |S| we see that when the maximum mass grows 
with A, the radius either decreases, or grows much more slowly. Overall, the stars 
that contain QM have radii between 8 and 12 km. 

As one would expect from Fig. |3] the heaviest pure NM stars (dotted lines) have 
larger radii >11 km, since the QM stars replace the NM stars at low radius (high 
central pressure). 


4 Conclusions 


We have seen that color superconductivity has a considerable effect on the equation 
of state for quark matter. We can see from Eq. (|2.14|) and Fig.|] (see also Ref. [jlT| ) 


that at values of the bag constant that would normally preclude compact stars from 
containing any QM, a large enough color-superconducting gap A can cancel out part 
of the bag constant, allowing a stable hybrid star to occur. In that sense, turning 
up A has a similar effect to turning down m s or turning down B , but there are 
differences (see below). 

The important question for the interpretation of future observations relates to 
the maximum mass of stars with quark matter in their core but with a surface made 
of nuclear matter, as observed recently [Q (see detailed discussion below). Our 
overall conclusion is that, in the range of bag constants and strange quark masses 
that we have studied, such stars have a maximum mass around 1.6 M 0 . This 
conclusion is based on our results presented in Fig. || The top panels show that 
color superconductivity allows QM-containing stars with a NM surface to occur at 
B 1//4 = 185 MeV (B = 153 MeV/frn 3 ), with masses up to about 1.5 M 0 . The lower 
panels show that at R 1 / 4 = 165 MeV (B = 96 MeV/frn 3 ) color superconductivity 
allows QM+NM stars to exist at high m s , and at low m s it boosts their mass: at 
m s = 150 MeV a nuclear-surface star of mass 1.6 M 0 is possible with A ~ 70 MeV, 
whereas the maximum mass would be 1.45 M 0 without color superconductivity. 

It should be borne in mind that the QM+NM stars with masses near 1.6 M 0 
are dominantly quark stars with a thin NM crust (see Fig. [|). This is what one 
would expect, given that they occur just to the left of the dots on the curves on 
Fig. | which mark the point at which the star becomes pure QM. The NM^QM 
transition in such stars occurs at very low pressure (< 1 MeV/fm 3 ) and at density 
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well below nuclear saturation density. 

It is also interesting to note (see Fig. |6j) that at low B and m s color 
superconductivity increases the maximum mass without appreciably changing the 
radius. This is different from the effect of changing the bag constant, which increases 
the mass and radius together. 

Finally, it is striking that even under the circumstances where color 
superconductivity has a noticeable effect on the maximum mass, it only does so 
for gap A > 50 MeV. Smaller gaps have little effect. 


We emphasize that we have had no difficulty in constructing stars that are 

El. 


compatible with the recent results of Cottam, Paerels, and Mendez 


who 


obtained a observational M-R curve by measuring the redshift of emission lines 
from highly ionized Oxygen and Iron in X-ray bursts from EXO0748-676. At radii 
of 8-12 km, their results suggest that the compact star mass is 1.2-1.8M 0 . Most 
of our stars fall in this range. Another constraint on the mass-radius of compact 
objects arise from recent observations of thermal radiation from RXJ185635-3754, 
an isolated neutron star. This combined with an accurate measurement of the 
distance to this object provide some information about its radius. However, since 
the spectrum is not quite black-body, the inferred radius depends on theoretical 
models employed to describe the objects atmosphere |26], For this reason, the 
constraint from RXJ185635-3754 is weak and not as compelling as those derived 
from EXO0748-676. A simple black-body fit yields radii that are small R ~ 5 
km |28|. In our study here we were unable to construct stars with such small radii. 
Atmosphere models which best fit the observed data favor a large radius. In a recent 
article, Walter and Lattimer find that these model studies indicate that R— 11.4 ±2 
km and M = 1.7 ± O.4M 0 |[7|. The superconducting quark stars constructed in 
this work are marginally compatible with these results. These larger values for the 
mass and radius require a small bag constant, a large A and small m s . 

We have shown that color superconductivity allows hybrid stars to have masses 
up to about 1.6 M 0 : does this mean that a definitive observation of a significantly 
heavier compact star would rule out quark matter? Our calculation could be 
improved in many ways (discussed below), so the upper limit we quote has 
theoretical errors that can only be roughly estimated. From Fig. [| we see that 
varying m s over the plausible range changes the maximum mass by less than 0.1 M 0 , 
and the mass rises from about 1.5 M 0 at R 1 / 4 = 185 MeV to about 1.6 M 0 
at R 1//4 = 165 MeV. Given this level of theoretical uncertainty, it seems that a 
definitive observation of a star with M > 1.8 M 0 would be difficult to explain in 
terms of hybrid QM+NM stars without invoking an even lower bag constant, with 
the danger that nuclear matter will be rendered unstable against two-flavor quark 
matter. 

There are many ways in which this line of inquiry could be pursued further. 
(1) We treated the quarks as free, with a color superconducting gap at their Fermi 
surface. We did not include perturbative corrections to the equation of state |29 
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in other words, we set the strong coupling constant a s = 0. It would be useful to 
perform such calculations for CFL quark matter, and see how robust our conclusions 
are against variation in a s . (2) We allowed two phases of quark matter, unpaired 
and CFL, and we did not include the two-flavor color superconducting “2SC” phase. 
Although this phase is generally unfavored f3(J, |3l|], it is just possible that there is 
a narrow range of m s in which it can occur. Another competitor is the crystalline 
phase [|32], |3 


It would be interesting to include these additional phases in our 


calculations. (3) We have used a general bag-model expression for the free energy of 
the quark matter, neglecting any possible density-dependence of the strange quark 
mass and color superconducting gap. We have also assumed that the bag constant 
takes the same value in all the quark matter phases, neglecting any differences of 
ground state energy between them. It would be interesting to use an NJL model to 
calculate the quark matter equation of state, and thereby include these phenomena. 
One could also use a bag constant with a phenomenological density dependence |33| . 

It is very encouraging that observations of the radii of compact stars are 
becoming more precise. There have been significant developments in the theory 
of dense quark matter in the last few years, and we look forward to seeing whether 
the observed properties of compact stars are compatible with (or even require) the 
presence of the exotic phases of quark matter that are being so widely discussed. 
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